The Onset of Phase Transitions in 
Condensed Matter and Relativistic QFT |^ 



R.J. Rivers, E. Kavoussanaki, 

Blackett Laboratory, Imperial College, London SW7 2BZ 

G. Karra. 

Salomon Smith Barney International Limited, London SWIW OSB 



Abstract 

Kibble and Zurek have provided a unifying causal picture for the appearance 
of topological defects like cosmic strings or vortices at the onset of phase transi- 
tions in relativistic QFT and condensed matter systems respectively. There is no 
direct experimental evidence in QFT, but in condensed matter the predictions are 
largely, but not wholly, supported in superfluid experiments on liquid helium. We 
provide an alternative picture for the initial appearance of strings/vortices that 
is commensurate with all the experimental evidence from condensed matter and 
consider some of its implications for QFT. 



^To be published in Condensed Matter Physics, 2000 
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1 The Problem 



One of the great successes in quantum field theory (QFT) in the last two decades has 
been the realisation of the unification of the disparate forces of quantum physics through 
spontaneous symmetry breaking. Even without a direct observation of the Higgs meson 
for electroweak unification, to date, supplementary evidence in support of the Standard 
Model of Glashow, Salam and Weinberg suggests that it is only a matter of time before 
Higgs particles will be found. 

Analogy with condensed matter physics would suggest that this and other symmetries 
were not always broken but that, in the very early, and very hot universe, they were 
restored. With several levels of unification anticipated we expect several changes of phase 
to have occurred, sequential partial orderings of an initially disordered state. Although 
the effects of most of the transitions can only be inferred indirectly, the last change 
of phase relevant to particle physics, the hadronisation of the quark-gluon plasma, is 
accessible at the heavy- ion collider facility at Brookhaven (RHIC), just about to begin 
to produce data. Although this transition is too difficult to address in this article, it 
provides a huge incentive to understand the way in which phase transitions occur in 
QFT. 

Consider a situation in which the symmetry group of the theory is broken, on cooling 
through the critical temperature Tc, by the degenerate groundstate manifold of the order- 
parameter fields. Using the tools of equilibrium thermal field theory (TFT) we can 
determine the nature of the transition. At late times after the transition the fields are 
ordered on large scales, in that they adopt a single value from this degenerate set over a 
large spatial region. 

To understand how this is achieved requires that we go beyond equilibrium TFT. In 
practice, we often know remarkably little about the dynamics of thermal systems. In 
particular, we shall restrict our discussion to the onset of such phase transitions, the very 
early times after the implementation of a transition when the scalar fields are only just 
beginning to become ordered. 

A simple question to ask is the following: In principle, the field correlation length 
diverges at a continuous transition. In practice, it does not. What happens? This is 
relevant for transitions that leave topological defects like walls, monopoles, vortices, or 
textures in their wake since we might expect 'defects' to be just that, entities whose sepa- 
ration is characterised by the correlation length. If this were simply so, an understanding 
of correlation lengths would lead to a prediction for defect densities. Conversely, a mea- 
surement of defect densities would be a measurement of correlation lengths. Vortices, 
in particular (cosmic strings), can be important for structure formation in the universe. 
Because of their implications for astrophysics, estimates of early field ordering in their 
presence have been made by KibbleJI], ||, using simple causal arguments. 

There are great difficulties in converting such predictions for the early universe into 
experimental observations. Zurek suggested that similar arguments were applicable 
to condensed matter systems for which direct experiments on defect densities, in par- 
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ticular vortex densities, could be performed. This has lead to a burst of activity from 
condensed matter experimentalists developing experiments to test these predictions. To 
date, agreement with the Zurek predictions for a variety of condensed matter systems is 
good, but not total. 

The problems are twofold, in that we need a better understanding of the predictions 
and of the experiments. In this article we shall only consider the first. Many of our 
observations have been discussed in more detail elsewhere, and our aim here is to give an 
overview. However, there is so much work in progress in this area that we do not claim 
to be exhaustive. 

2 The Predictions 

The predictions of Kibble and Zurek for the way in which the fields freeze in at the onset 
of a transition are simply expressed in terms of causality. As such, they are generic. 
However, to quantify them it is convenient to consider simple field theories in which they 
can be realised. 

In principle, the simplest transitions are temperature quenches. In practice, the 
circumstances in which transitions occur are poorly understood. In the early universe 
we do not necessarily begin from a thermal distribution. More awkwardly, in condensed 
matter physics (and heavy-ion collisions) we do not have a spatially uniform temperature. 
However, for the sake of argument we assume a uniform temperature, and leave more 
realistic scenarios until later. 

There is a further problem. In condensed matter systems the nature of the fields in the 
Landau-Ginzburg effective theory is prescribed, and their number is small. In the early 
universe we have little constraint on our model making beyond a feeling that, despite the 
ubiquity of gauge fields at later times, at early times scalar fields are important. 

2.1 The models 

As we shall see later, most experiments involve vortices and, with this in mind, we adopt 
the simplest theory, that of a single complex scalar field in three spatial dimensions 
with a wine-bottle potential. The transition is continuous^]. 

That is, we assume that the qualitative dynamics are conditioned by the field's equi- 
librium free energy, of the form 

F(T)= Jd'x ^\V<P\' + m\T)\<pf + X\<P\'y (2.1) 

Intentionally, with the early universe in mind we have written F{T) of Eq. |2.1| in the form 
appropriate for a relativistic quantum field. The coefficient m^(T) has the interactions 

^There is no evidence for strong first order, or discontinuous, transitions in QFT. The mechanism for 
such transitions, bubble nucleation, would lead to a very different onset than that considered here. 



3 



with the real particles of the heat bath taken into account, and vanishes at T = Tc- 
We are considering the field as an open system, in which changes in the external 
environment lead to changes in the parameter m^(T), which takes the value — when 
T = 0. For relativistic QFT the change in temperature that leads to the change in the 
sign of rn? is most simply understood as a consequence of the system expanding. Thus, 
in the early universe, once thermalisation is possible, a weakly interacting relativistic 
plasma at temperature T ^ M has an entropy density s oc T^. As long as thermal 
equilibrium can be maintained, constant entropy S per comoving volume, S oc sa{t)^, 
gives T oc a{t)~^ and falling, for increasing scale factor a{t). 
Specifically, at time t the temperature T{t) satisfies 

where mp is the Planck mass, and N* is the number of effective field degrees of free- 
dom. Models that attempt to take infiation into account, however, lead to 'preheating' 
that is not Boltzmannian[^]. Nonetheless, even in such cases it is possible to isolate an 
effective temperature for long-wavelength modes. This is all that is necessary, but is too 
sophisticated for the simple scenarios that we shall present here. Even the inclusion of 
an FRW metric would complicate the issue at this stage, and we assume flat space-time, 
with decreasing temperature. Prior to the transition, we assume a uniform temperature 
T > Tc, for which m^(T) > 0. After the transition, m^(0) = — < enforces the U{1) 
symmetry-breaking, with fleld expectation values (0) = ±?7, t]"^ = M^/ A. The Compton 
wavelength = is the natural distance scale. For the sake of argument we assume 
mean-fleld behaviour m?{T) = M^(T/Tc — 1), whereby the equilibrium correlation length 
UT) = |m(T)|-i = eo(T/T,-l)-V2. 

Equally well, after rescaling, F could be the Ginzburg-Landau free energy 

F(T) = Jd^x (^|!-|V0p + a(T)|0p + /3|0r) (2.3) 

for a non- relativistic condensed matter fleld, in which the chemical potential tt(T) = 
ao{T/Tc — 1) vanishes again at the critical temperature Tc. In this case we envisage 
the change in a{T) as a consequence of an external cooling of the system or a change 
in the pressure of the system that leads to a change in T^. In either case, we again 
assume circumstances in which, in a flnite time, T/Tc varies from greater than unity to 
less. The fundamental length scale is given from Eq.^]3| as = h /2mao whereby 
ieq{T) = ^q{T /Tc — as before. The Gross-Pitaevski theoryQ suggests a natural 

time-scale tq = Ti/aQ. When, later, we adopt the time- dependent Landau- Ginzburg 
(TDLG) theory we flnd this still to be true, empirically, at order-of-magnitude level, and 
we keep it. 

The minima of the flnal potential now constitute the circle = ?7e*", where rf = M^/ A 
for Eq. |2.1| , and rj"^ = ao/jS for Eq. |2.3| . When the transition begins begins to fall 
into the valley of the potential, choosing a random phase. Although this randomly 
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chosen phase will vary from point to point we expect domains across which the phase is 
roughly constant. How this collapse takes place determines the size of the first identifiable 
domains. It was suggested by Kibble and Zurek that this size is essentially the equilibrium 
field correlation length at some appropriate temperature close to the transition. Two 
very different mechanisms have been proposed for estimating this size. 

2.2 A first guess: Thermal activation 

In early work on transitions it was assumed [|l| that initial domain size was fixed in the 
Ginzburg regime, identifiable from equilibrium theory. By this we mean the following. 
Suppose the temperature T(t) varies sufficiently slowly with time t that it makes sense to 
replace V{(f),T) by ¥{(/), T{t)). Well away from the transition this is justified, but close 
to the transition it is not. Once we are below T^, and the central hump in \^(0,T(t)) is 
forming, the Ginzburg temperature Tq < Tc signals the temperature above which there is 
a significant probability for thermal fluctuations over the central hump on the scale of the 
correlation length at that temperature. Most simply, it is determined by the condition 

AV{TG)ee,iTG) ^ ksTa (2.4) 

where AV{T) is the difference between the central maximum and the minima of V^(0, T). 
That is, 

U^g)=0(^^-^^), (2.5) 

Whereas, above Tq there will be a population of 'domains', fluctuating in and out of 
existence, at temperatures below Tq fluctuations from one minimum to the other become 
increasingly unlikely. For the relativistic theory of Eq. p.l| , in units in which fc^ is unity, 
we find |1 — Tc/Td = 0(A), whereby 

On the other hand, for non-relativistic condensed matter, we find 

1 - Tg/T, = {PkBTjalilf (2.7) 

It was originally suggested by Kibble that we identify ^eqiTc) with the scale at 
which stable domains begin to form. We shall show this to be incorrect, for quenches 
that are not too slow. However, we will find that strong thermal fluctuations do have a 
role, particularly in ^ife, for which the whole experiments take place within the Ginzburg 
regime. This is an issue that requires more than equilibrium physics. The most simple 
dynamical arguments invoke causality. 
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2.3 Causality in QFT and condensed matter 



The first application of causality was again due to Kibble 0. Whatever the details of the 
microscopic physics, causality sets an upper limit over which the field can be correlated 
in the causal horizon with diameter r{t) ~ 2t. In the vicinity of the transition, whether 
at T = Tg or not, this gives a correlation length ^(t) < r{t) or, from Eq. p.2| 



where N* the number of field degrees of freedom at T^. This was used to bound monopole 
density in the early universe. Although the lack of correlation in the field phase in 
different causal horizons demonstrates the existence of a field structure that will naturally 
lead to defects, beyond that it is not a helpful guide. 

If this invocation of causality, or the Ginzburg criteria, attempt to set scales once 
the critical temperature has been passed, other causal arguments attempt to set scales 
before it is reached. Again suppose that T/Tc varies in time as a consequence of the 
change in the environment. We have seen that ^eq{T{t)), obtained by inserting the time- 
dependence of T/Tc explicitly, diverges at T(t) = Tc, which we suppose happens at t = 0. 
This cannot be the case for the true correlation length ^(t), which can only grow so far 
in a finite time. Initially, for t < 0, when we are far from the transition, we again assume 
effective equilibrium, and the field correlation length C,(t) tracks ^eq{T{t)) approximately. 
However, as we get closer to the transition, ^eq{T{t)) begins to increase arbitrarily fast 
and the adiabatic approximation breaks down. This is the largest value that attains. 
This largest value prior to the transition corresponds to the smallest value of the field 
correlation length after the transition. We shall argue later that the argument is too 
simple but, nonetheless it is a plausible starting point. The question then becomes: how 
large does the field correlation get in practice? Whether in QFT or condensed matter 
physics there is a maximum speed at which ^ (t) can grow on purely causal grounds. As a 
crude upper bound, in QFT the true correlation length ^(t) fails to keep up with ^eq{T(t)) 
by the time —i at which ^^q is growing at c=l, the speed of light, d^eq{T{—t))/dt = 1. 
It was suggested, again by Kibble [0], that once we have reached this time ^{t) freezes 
in, remaining approximately constant until the time t ~ +t after the transition when 
it once again becomes comparable to the now decreasing value of ^eg. The correlation 
length ^ = ^eqij) = ^eq{—t) is argued to provide the scale for the minimum domain size 
after the transition. 

Specifically, if we assume a time-dependence m^(t) = —M'^t/tq in the vicinity of 
t = 0, when the transition begins to be effected, then the causality condition gives 




(2.8) 



tc = tU\2M)-'^/\ As a result. 



M^eq{i) 



(2.9) 



which we write as 




1/3 



(2.10) 
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where tq = = are the natural time and distance scales. In contrast to Eq. p.5| , 
Eg. [2.101 depends explicitly on the quench rate, as we would expect. For tq ^ tq the field 
is correlated on a scale of many Compton wavelengths. 

This approach of Kibble was one of the motivations for a similar analysis by Zurekj^ of 
transitions with scalar order parameters in condensed matter. Explicitly, in the Ginzburg- 
Landau free energy Eq. |2.3| , <y{T) also vanishes at the critical temperature T^. The only 
difference is that, in the causal argument, the speed of light should be replaced by the 
speed of (second) sound, with different critical index. 

Explicitly, let us again assume the mean-field result «(T) = ao^iT), where e = 
{T/Tc — 1), remains valid as T/Tc varies with time t as a{t) = a{T{t)) = —aot/Tg in the 
vicinity of T^. It follows that the equilibrium correlation length ^eq{t) and the relaxation 
time r(t) diverge when t vanishes as 



t 



-1/2 



r{t) = To 



t 



(2.11) 



in terms of and Tq as given earlier. 

The speed of sound is c(t) = ^f,g(t)/T(t), slowing down as we approach the transition 
as The causal counterpart to d^eq{t)/dt = 1 for the relativistic field is d^eq{t)/dt = 

c{t). This is satisfied at t = — t, where t = ^tqTq, with corresponding correlation length 



e = Ut) = U~t) = ^o[-^) . (2.12) 

(cf. Eq. p.lOl ). A variant of this argument]^ that gives essentially the same results 
is obtained by comparing the quench rate directly to the relaxation rate of the field 
fiuctuations. We stress that, yet again, the assumption is that the length scale that 
determines the initial correlation length of the field freezes in before the transition begins. 
Whatever, the field is already correlated on a scale of many Compton wavelengths when 
it begins to unfreeze. 



2.4 Experimental predictions 



The end result of the simple causality arguments is that, both for QFT and condensed 
matter, when the field begins to order itself its correlation length has the form 



(2.13) 



for appropriate 7. In fact, the powers of Eq. p.lO] and Eq. |2.12| are mean- field results, 
changed on implementing the renormalisation group. Whereas, for ^He, the critical 
behaviour of Eq. |2.12| survives, for ^He, 7 = 1/3. On the other hand, for QFT the 
relevant parameter is ^T^, the ratio of the maximum correlation length to the thermal 
length = T~^. In equilibrium theory, when the temperature is high enough the theory 
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is essentially three-dimensional, with critical index 7 different from its /onr-dimensional 
mean-field value. At our level of discussion, it is sufficient to keep (/owr-dimensional) 
mean-field values. 

As we said earlier, when the transition begins begins to fall into the valley of the 
potential, choosing a random phase. This randomly chosen phase will vary from point to 
point leading to approximate domains across which the phase is roughly constant. Such 
domains will meet at defects, in this case vortices, tubes of 'false' vacuum 0^0, around 
which the field phase changes by ±2-7?. In an early universe context these are 'cosmic 
strings', but we shall not consider their properties here. 

Correlation lengths in the early universe are not amenable to direct observation. 
Kibble made a second assumption, that the correlation length Eq. p.lOj also sets the scale 
for the typical minimum intervortex distance at the time that vortices are produced. 

That is, the initial vortex density riaef is 

for 7 = 1/3, where / = 0(1) estimates the fraction of defects per 'domain'. Equivalently, 
the length of vortices in a box volume v is 0{ndefv). 

What is striking and suspect about these predictions is that they are universal. They 
do not use any information about the strength (A or (3) of the interactions, and hence the 
magnitude of the order parameter after the transition or, in consequence, the existence 
of the Ginzburg regime. 

Since ^0 also measures cold vortex thickness, tq ^ tq corresponds to a measurably 
large number of widely separated vortices. For the early universe Kibble deduced 



^=i—=^ , (2.15) 

where, as before, mp is the Planck mass, and A^* the number of field degrees of freedom 
at Tc0. If it could be argued that this initial network behaves classically then, thereafter, 
the density will reduce due to the collapse of small loops, intersections chopping off loops 
which in turn collapse, and vortex straightening so as to reduce the gradient energy of 
the field. 

However, even if cosmic strings were produced in so simple a way in the very early 
universe it is still not possible to compare the density Eq. |2.14| with experiment. What 
is more amenable to experiment, in principle, is the length distributions of string net- 
works, and their ability to show scaling behaviour. This only impinges indirectly on the 
correlation length of the field, but would have to be commensurate with any density 
calculations. 



■^We note that, despite their different origins, the predictions 2.6, 2. 8, 2. f 5 may not differ so hugely 
from one another for very smaU couphng. 
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It was Zurek[0] who first suggested that, if this relationship Eq. |2.14| between defect 
density and correlation length were true, it could be tested directly in condensed matter 
systems, particularly in liquid helium. 



3 Experiments in Condensed Matter 
3.1 Vortices in superfiuid helium 

Vortex lines in both superfiuid '^He and ^He are analogues of global cosmic strings. A 
crude but effective model is to treat the system as composed of two fluids, the normal fluid 
and the superfiuid, which has zero viscosity. In ^ife the bose superfiuid is characterised 
by a complex field 0, whose squared modulus |0p is the superfiuid density. The superfiuid 
fraction is unity at absolute zero, falling to zero as the temperature rises to the lambda 
point at 2.17K. The Landau-Ginzburg theory for ^ife has, as its free energy F{T) of 
Eq.p. 

The situation is more complicated, but more interesting, for ^ife, which becomes su- 
perfiuid at the much lower temperature of 2 mK. The reason is that the ^He is a fermion. 
Thus the mechanism for superfluidity is very different from that of '^He. Somewhat as in 
a BCS superconductor, these fermions form the counterpart to Cooper pairs. However, 
whereas the (electron) Cooper pairs in a superconductor form a state, the ^He pairs 
form a state. The order parameter Aai is a complex 3x3 matrix Aai- There are 
two distinct superfiuid phases, depending on how the S0{3) x 5*0(3) x f/(l) symmetry 
is broken. If the normal fluid is cooled at low pressures, it makes a transition to the 
^He — B phase, in which A^i takes the form A^i = Rai{uj)e^'^, where i? is a real rotation 
matrix, corresponding to a rotation through an arbitrary (^0 0. 

The Landau-Ginzburg free energy is, necessarily, more complicated0, but the effec- 
tive potential V{Aai,T) has the diagonal form V{A,T) = a(T)\Aai\'^ + 0{A'^) for small 
fluctuations, and this is all that we need for the production of vortices at very early 
times. Beyond that it can be mimicked by Eq.p73| for our purposes (e.g. see 0). Thus 
the Zurek analysis leads to the prediction Eq. p.l4| , as before, for appropriate 7. However, 
for ^He the mean- field approximation is good and the mean-field critical index 7 = 1/4 is 
not renormalised, whereas for ^ife a better value is 7 = 1/3, as for the naive relativistic 
theory. 



3.2 Experiments in ^He. 



Although ^He is more complicated to work with, the experiments to check Eq. |2.14| are 



cleaner in that, because the nucleus has spin 1/2, even individual vortices can be detected 
by magnetic resonance. Further, because vortex width is many atomic spacings the 
Landau-Ginzburg theory is reliable. 



*At large distance scales there is a complication in that the small spin-orbit coupling becomes im- 
portant, to fix oj at arccos(— 1/4), but this will not concern us here. 
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So far, experiments have been of two types. In the Helsinki experiment]^ superfluid 
^He — 5 in a rotating cryostat is bombarded by slow neutrons. Each neutron entering the 
chamber releases 760 keV, via the reaction n+^He He + 7Q0keV. The energy goes 

into the kinetic energy of the proton and triton, and is dissipated by ionisation, heating 
a region of the sample above its transition temperature. The heated region then cools 
back through the transition temperature, creating vortices. Vortices above a critical size 
(dependent on the angular velocity of the cryostat) grow and migrate to the centre of the 
apparatus, where they are counted by an NMR absorption measurement. Suffice to say 
that the quench is very fast, with tq/tq = O(IO^). Agreement with Eq. |2.14| and Eg . p. 12 



is very good, at the level of less than an order of magnitude. This is even though it is 
now argued IQ that the Helsinki experiment should not show agreement because of the 
geometry of the heating event. 



The second type of experiment has been performed at Grenoble and Lancaster 
Rather than count individual vortices, the experiment detects the total energy going 
into vortex formation. As before, ^He is irradiated by neutrons. After each absorption 
the energy released in the form of quasiparticles is measured, and found to be less than 
the total 760 keV. This missing energy is assumed to have been expended on vortex 
production. Again, agreement with Zurek's prediction Eg. 2. 14 and Eg. |2.12 is good. 



3.3 Experiments in '^He. 

The experiments in ^ife, conducted at Lancaster, follow Zurek's original suggestion. The 
idea is to expand a sample of normal fluid helium, in a container with bellows, so that 
it becomes superfluid at essentially constant temperature. That is, we change 1 — T/Tc 
from negative to positive by reducing the pressure, thereby increasing Tc. As the system 
goes into the superfluid phase a tangle of vortices is formed, because of the random 
distribution of field phases. The vortices are detected by measuring the attenuation of 
second sound within the bellows. Second sound scatters off vortices, and its attenuation 
gives a good measure of vortex density. A mechanical quench is slow, with tq some tens of 
milliseconds, and tq/tq = O(10^°). Two experiments have been performed [pJ], |12|. In the 



first fair agreement was found with the prediction Eg. |2.14| , although it was not possible 
to vary tq. However, there were potential problems with hydrodynamic effects at the 
bellows, and at the capillary with which the bellows were filled. A second experiment, 
designed to minimise these and other problems has failed to see any vortices whatsoever. 

Some care is needed. Not only is the Landau-Ginzburg effective theory more suspect 
for ^He, but its Ginzburg regime is so wide, at 0{1K) that the transition takes place 
entirely within it. 

3.4 Quenching in an annulus 

Experiments for vortex densities are problematical in that a closely bound vortex-antivortex 
pair gives a count of 2 right up to annihilation, whereas their topological charge remains 
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zero. It is possible to devise experiments that count topological charge. Consider a closed 
path in the bulk superfiuid with circumference C ^ .^(t). Naively, the number of 'regions' 
through which this path passes in which the phase is correlated is TV = 0{C/^{t)). As- 
suming an independent choice of phase in each 'region', the r.m.s phase difference along 
the path is 



A9c-VAf = O{^C/0. (3.16) 

If we now consider a quench in an annular container of similar circumference C of 
superfiuid '^He and radius / ^ C, Zurek suggested^ that the phase locked in is also 
given by Eq. |3.16| , with ^ of Eq. |2.12| . Since the phase gradient is directly proportional to 



the superflow velocity we expect a flow after the quench with r.m.s velocity 

Av = o(-.^). (3.17) 



provided / = O(^). Although in bulk fluid this superflow will disperse, if it is constrained 
to a narrow annulus it should persist, and although not large is measurable, in principle. 
Specifically, in the units of Zurek [|| 



Av = 0((rQ[H)-'^/7\/cM) (3.18) 

where tq is, typically, tens of milliseconds and C in centimetres, z/ = 1/2 is the mean- 
field critical exponent above. In principle u should be renormalised to u = 2/3, but the 
difference to Av is negligable. With such a small index the result is almost independent 
of quench rate. 

In practice there are difficulties in performing annular measurements in ^ife. A 
similar, but easier, experiment can be performed on annular Type-II superconductors, 
on cooling through their critical temperatures. The relevant free energy is|p the extension 
of F of Eq.p|, 

F(T)= Jd'x ^±.\-^h\/<P-^A\' + aiTM' + ^P\M'^+^. (3.19) 

A is the vector potential in the Coulomb gauge, and B = V A A. 

If we can initially ignore the effects of the gauge field as a temperature quench is 
imposed through a change in a{T), as before, the result Eq. |3.16| persists. For perimeter 



C the variance in the number of flux quanta produced spontaneously is 

A^c=^A«c-^/f. (3.20) 

In fact, it is more convenient to quench annular Josephson Junctions, in which two 
identical rings of superconductor are held apart by an oxide layer through which Cooper 
pairs can tunnel. If 6i and 62 are the phases of on the upper and lower rings then, once 
the transition has taken place the tunneling current has the form 

J = JcSin{9i - 62) (3.21) 
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and the theory is described by a dissipative Sine-Gordon equation. 

The kinks of this equation are the 'fluxons' of the Josephson Junction and are easy 



to observe experimentally! 13]. The variance in fluxon number at their formation is 



ANc = ^A{ei-e2). (3.22) 
Zn 

Since 61,62 are independent as the transition begins, we would expect, from Zurek's 
analysis, that 



whence 




(3.23) 



(3.24) 



Although some caution is required in the interpretation of the experiments, which were 
not devised with this prediction in mind, it seems to be supported by experiment, for 
which, with ^ = O(10~^)mm and C = 0.5xl0~^mm, say, we would expect AA^^ = 0(1). 

There is certainly no agreement, in this or any other experiment in which defects are 
seen, with the thermal fluctuation density that would be based on Eq.|! 



4 The Kibble-Zurek picture for the freezing in of ^ 
is correct. 

Since all equations of motion have causality built into them we should be able to confirm 
the first predictions Eq. 2.1(]| and 2.12| of Kibble and Zurek explicitly, as we shall now see. 



4.1 Condensed matter: the TDLG equation 

We assume that, for the condensed matter systems of interest to us, the dynamics of 
the transition can be derived from the explicitly time- dependent Landau-Ginzburg free 
energy 

m = J d'x {^WaY + ami + ■ (4.25) 

in which we substitute T{t) for T directly in (gj). In (lOSi ) = (0i + #2)7^2 (a = 1, 2) 
is the complex order-parameter field, whose magnitude determines the superfiuid density. 
As before, in a mean field approximation, the chemical potential a(T) takes the form 
a(T) = aQe{t), where e = {T/T^ — 1). In a quench in which or T changes it is 
convenient to shift the origin in time, to write e as 

e{t) = eo - —6{t) (4.26) 



5 Corresponding to tq = 0(l)s, $0 = O{10^)A. 



12 



for — oo < t < rQ(l + eo), after which e(t) = — 1. eo = 1 — Tq/Tc measures the original 
temperature Tq and tq defines the quench rate. The quench begins at time t = but the 
transition from the normal to the superfluid phase only begins at time to = ^o^q- When 
it is convenient to measure time from the onset of the transition we use the notation 
At = t- to. 



Motivated by Zurek's later numerical ||14]| simulations, we adopt the time-dependent 
Landau- Ginzburg (TDLG) equation for F, 

(4.27, 



r dt S^a 

where rja is Gaussian thermal noise, satisfying 

(r/,(x, t)r/,(y, t')) = 25„,T(t)r5(x - y)S{t - t'). (4.28) 

This is a crude approximation for ^He, and a simplified form of a realistic description of 
^He but it is not a useful description of QFT, as it stands. 

It is relatively simple to determine the validity of Zurek's argument since it assumes 
that freezing in of field fluctuations occurs just before symmetry breaking begins. At that 
time the effective potential V{(f),T) is still roughly quadratic and we can see later that, 
for the relevant time-interval —t < At < i the self-interaction term can be neglected 
(/3 = 0). 

In space, time and temperature units in which C,o = '^o = = 1, Eq. [4.27| then 
becomes 

0,(x, t) = -[-V^ + e(t)]0„(x, t) + r/,(x, t). (4.29) 

where f] is the renormalised noise. The solution of the 'free'-field linear equation is 
straightforward, giving a Gaussian equal-time correlation function 

(0a(r,t)0,(O,t)) = 6atG{r,t) = J fke'^'-'Pik^t). (4.30) 

in which the power spectrum P{k, t) has a representation in terms of the Schwinger 
proper-time r as 

P(k,t) = r dTfit-T/2)e-'^^\-^^o^''^'-''''\ (4.31) 



where T is the renormalised temperature. In turn, this gives p5 



oo 



/ 1 \ 3/2 r 

G{r,t)= I dTf{t-T/2) (^j e"'' /^"e-Jo (4.32) 







For constant e, as happens at early times, on rescaling in Eq. [4.32| we recover the usual 
Yukawa correlator 

G(r,t) = |Le-'^/?o, (4.33) 
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where Tq = Tc(l + eo) is the initial temperature. For e(t) of Eq. [4.26| a saddle-point 
calculation shows that, although we recover Eq. [4.33| for r/^> {^/^qY later times, the 
correlation function is then dominated by its smaller-r behaviour. At time to = eoTo, 
when the transition begins, a saddle-point calculation shows that, provided the quench 
is not too fast. 



G(r,^o)^^e-'^^^/^^^^^ (4.34) 

where a = 0(1), confirming Zurek's result. 

Zurek's prediction is robust, since further calculation shows that ^(t) does not vary 
strongly in the interval —i<At< i, where At = t — to- 



4.2 QFT: Closed time-path ensemble averaging 

For QFT the situation is rather different. In the previous section, instead of working with 
the TDLG equation, we could have worked with the equivalent Fokker-Planck equation 
for the probability pt[^] that, at time t > 0, the measurement of will give the function 
$(x). When solving the dynamical equations for a hot quantum field it is convenient to 
work with probabilities from the start. 

Take t = as our starting time for the evolution of the complex field = (0i + 
%(\)2)j^2. Suppose that, at this time, the system is in a pure state, in which the mea- 
surement of would give $o(x). That is:- 

0(t = 0,x)|$o,t = 0) = <l>o|$o,t = 0). (4.35) 

The probability that, at time tj > 0, the measurement of will give the value $ is 
j9t[$] = I^E'oP, where \E'o is the state-functional with the specified initial condition. As a 
path integral 

^0= / l)0exp zSt[0] , (4.36) 

where 5*4 [0] is the (time-dependent) action that describes how the field is driven by 
the environment and spatial and field labels have been suppressed (e.g. P0 = P0iP02). 
Specifically, for t > the action for the field is taken to be 

SM\ = I dx(^^d,cP,d>^<P^ - ^m'ml - Ixi^Pir). (4.37) 

where m(t) describes the evolution of the action under external influences, to which the 
field responds. 

It follows that can be written in the closed time-path form 

= P0+P0_ exp{t(s,[<P+] - 5,[0_]) j, (4.38) 

J^±(0)=<I.o I V /J 

where 1^0^ = 'P0±,i^</'±,2- Instead of separately integrating 0-|- along the time paths 
< t < tj, the integral can be interpreted as time-ordering of a field along the closed 
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path C+ ® C_ where = 0+ on C+ and = 0_ on C_. When we extend the contour 
from tf to t — oo either 0+ or 0_ is an equally good candidate for the physical field, but 
we choose 0+. 







a 



+ 



Figure 1: The closed timepath contour C+ ® C_. 



The choice of a pure state at time t = is too simple to be of any use. As we 
said earlier, we assume that $ is Boltzmann distributed at time t = at an effective 
temperature of Tq — according to the Hamiltonian H[^] corresponding to the action 
S'[0], in which is taken to be periodic in imaginary time with period Pq. We now have 
the explicit form for pt [$] , 



written as the time ordering of a single field along the contour C — C+ ® C_ ® C3, 

extended to include a third imaginary leg, where takes the values 0+, 0_ and 03 on 
C+, C_ and C3 respectively, for which Sc is >S'[0+], 5'[0_] and 5*0 [03]. 




(4.39) 



to 



a 



+ 



C- 



3?e(r) 



C3 



to - if^o 



Figure 2: A third imaginary leg 
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To demonstrate how we can average without having to calculate pt[$] explicitly we 
see that Gab{\^ - x'|; t) = ($a(x)$f,(x'))t is given by 

Gati\^-^'\;t) = (0,(x,t)0,(x',t)), (4.40) 

the equal-time thermal Wightman function with the given thermal boundary conditions. 
Because of the non-equilibrium time evolution there is no time translation invariance in 
the double time label. 



4.3 QFT: the free roll 



Fortunately, as for the condensed matter case, the interval —i < At < i occurs in the 
linear regime, when the self-interactions are unimportant. The relevant equation for 
constructing the correlation functions of this one-field system is now the second-order 
equation 

for F of Eq. p.lj . This is solvable in terms of the mode functions Xk{^)i identical for 
a = 1,2, satisfying 

4 + + 



(4.42) 



subject to Xfc(^) 



at t < 0, for incident frequency uOi, 



+ eoM2, for 



m?{t) = e{t)M'^, where e(t) is parameterised as for the TDLG equation above. This 
corresponds to a temperature quench from an initial state of thermal equilibrium at 
temperature Tq > Tc, where {Tq/Tc — 1) = eo- There is no reason to take eo small. The 
diagonal correlation function G{r,t) of Eq. [1.30| is given as the equal-time propagator 



G{r,t) 



27r2 



fke'^-xtm,:it)G{k) 



dke'-^xtm-,it)Gik), 



(4.43) 



kr 



where G{k) = coth.{ujin{k) /2To)/2u!in{k) encodes the initial conditions. 

An exact solution can be given |T6[] in terms of Airy functions. Dimensional analysis 
shows that, on ignoring the k-dependence of G{k), appropriate for large r (or small k), 
ieqit) of Eq. |2.1CI| again sets the scale of the equal-time correlation function. Specifically, 



/sm 
dK — 
K 



sin /t(r/^) 



F{k), 



(4.44) 



where -F(O) = 1 and F{k) ~ k for large k. Kibble's insight is correct, at least for this 
case of a single (uncoupled) field. 
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5 However, Defect Densities do not Determine ^ Di- 
rectly (and Vice- Versa) 

We have seen that there is no reason to disbeheve the causal arguments of Kibble for 
QFT and Zurek for condensed matter as to the correlation length ^ at the onset of 
the transition. The excellent agreement with the ^He experiments suggests that, for 
condensed matter, this length does, indeed, characterise vortex separation at the time 
when the defects form. 

However, if we take the Lancaster experiment at face value, this cannot always be the 
case. This is not surprising. Numerical simulations of TDLG equations, by Zurek himself, 
and the related Langevin equations for non-equilibrium statistical fields show that, at 
early times after a transition, the field fluctuations remain on very small scales. In reality, 
different frequency modes freeze at different times, and the causal argument only applies 
strictly to the long wavelength modes. Of course, field ordering is controlled by the long 
wavelength modes, but we would like a more complete description of the freezing in of 
the field, if we are to make the second step of relating this to defect structure in the 
field . That is, although the field is correlated on long scales ^, this does not preclude 
structure on very much shorter scales. Moreover, it is this structure, that we might wish 
to think of as proto-defects, that will evolve into the defects, in this case vortices, once 
the transition is complete. Inspection of the field at early times shows no obvious relation 
between the separation of these proto-defects and ^. The question then becomes one of 
why the ^He experiments should be in agreement, rather than why the ^ife experiments 
are not. 



5.1 Classical vortices in condensed matter and QFT 

The 0(2) string is, classically, a tube of false vacuum of radius 0{M~^) whose core is 
characterised by a line of string zeroes. Its winding number n G Z is the change in phase 
of the complex field around the core (in units of 27r). We shall only consider strings with 
|n| = 1 since all others can be considered as multiple zeroes. 

It would be foolish to estimate the probability of finding vortices directly from J9t[<l']. 
A starting-point for counting vortices in superfiuids is to count line zeroes. Not all line 
zeroes are candidates for vortices since zeroes occur on all scales. However, a starting- 
point for counting vortices in superfiuids is to count line zeroes of an appropriately coarse- 
grained field, in which structure on a scale smaller than ^o? the classical vortex size, is not 
present [ffTf. This is also the unstated basis of the numerous numerical simulations [ITS] of 



cosmic string networks built from Gaussian fluctuations (but see |]T9[)- 

Even then, there are several prerequisites before line zeroes can be identified with 
vortex cores, and n^eroit) with ride fit). 

• The field, on average must have achieved its symmetry-broken ground-state equi- 
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librium value 

= ao/P orMVA, (5.45) 

non-perturbatively large (in /3) This, in itself, is sufficient to show that the causal 
time t is not the time to begin looking for defects, since is small at this time. 

• Only when drizero/dl is small in comparison to Uzero/^ at / = "will the line-zeroes 
have the small-scale non-fractal nature of classical defects, although defects may 
behave like random walks on larger scales. As the power in the long wavelength 
modes increases the 'Bragg' peak develops in k'^G{k,t), moving in towards k = 0. 
This condition then becomes the condition that the peak dominates its tail. 

• The energy in field gradients should be commensurate with the energy in classical 
vortices with the same density as that of line zeroes. 

We stress that these are necessary, but not sufficient, conditions for classical vortices. 
In particular, only the full nonlinearity of the system can establish classical profiles. 
We will term such zeroes as satisfy these conditions, proto- vortices. In fact, most (but 
not all[^]) numerical lattice simulations cannot distinguish between proto- vortices and 



classical vortices. 

Whereas the above are equally true for condensed matter and QFT there are further 
complications peculiar to QFT. In particular, in QFT we need to consider the whole 
density matrix ($'|p(t)|$) rather than just the diagonal elements Pt[^] = ($|p(t)|$). 
Classicality is understood in terms of 'decoherence', manifest most simply by the ap- 
proximate diagonalisation of the reduced density matrix on coarse-graining. By this we 
mean the separation of the whole into the 'system', and its 'environment' whose degrees 
of freedom are integrated over, to give a reduced density matrix. The environment can 
be either other fields with which our scalar is interacting or even the short wavelength 
modes of the scalar field itself [^1], |2^. When interactions are taken into account this 
leads to quantum noise and dissipation. 

In the Gaussian approximations for QFT that we shall adopt here, with ($) = 0, 
integrating out short wavelengths with k > is just equivalent to a momentum cut- 
off at the same value. This gives neither noise nor dissipation and diagonalisation does 
not occur. Nonetheless, from our viewpoint of counting line-zeroes, fluctuations are still 
present when / = 0{M~^) that prevent us from identifying line- zeroes with proto- vortices 
easily. 

For all these caveats, there are other symptoms of classical behaviour in QFT once 
Gi{0;t) is non-perturbatively large. Instead of a field basis, we can work in a particle 
basis and measure the particle production as the transition proceeds. We shall see that, 
in the Gaussian approximation, n^ero, the density of line zeroes, is given in terms of the 
moments of Gi{r, t) . Whether we expand with respect to the original Fock vacuum or with 
respect to the adiabatic vacuum state, the presence of a non-perturbatively large peak 
in k'^G{k]t) at k = ko signals a non-perturbatively large occupation number Nkg oc 1/A 
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of particles at the same wavenumber fco[^. With rtzero of order /cq this shows that the 
long wavelength modes can now begin to be treated classically. 

From a slightly different viewpoint, the Wigner functional only peaks about the clas- 



sical phase-space trajectory once the power is non-perturbatively large [P4| , |25| . More 
crudely, the diagonal density matrix elements are only then significantly non-zero for 
non-perturbatively large field configurations oc A~^/^ like vortices. 

5.2 Line-zero density 

Suppose, at some time, that the field has line zeroes R„(s), where n = 1,2, .. labels the 
zero, and s measures the length along it. As a result the topological line density of zeroes 



p(r) can be defined p6|, [27[] by 



p(x) = Y.j ds^S'l^ - Rn(3)]. (5.46) 



In ( ^.46|) ds is the incremental length along the line of zeroes Rn(s) (n=l,2,.. .) and 



is a unit vector pointing in the direction which corresponds to positive winding number. 
It follows that, in terms of the zeroes of the field $(x), Pi(x) can be written as 

p,(x) = 52[$(x)]e,,,9,$i(x)9fc$2(x), (5.47) 

where 5^[$(x)] = 5[$i(x)](5[$2(x)]. The coefficient of the 5-function in ( p.47| ) is the 
Jacobian of the more complicated transformation from line zeroes to field zeroes. What 
we want is not this, but the total line density p(x), 

p,(x) = 52[$(x)]|e,,fca,<l'i(x)afc$2(x)|. (5.48) 

The vanishing field expectation value and the independence of the field and its deriva- 
tives 

(<I>„(x)) = = ($,(x)a,$;,(x)), (5.49) 

imply (pj(x)) =0 i.e. an equal likelihood of a string line- zero or an antistring line-zero 
passing through an infinitesimal area. However, 

n{t) = (A(x))i > (5.50) 

and measures the total line-zero density in the direction i, without regard to string 
orientation. The isotropy of the initial state guarantees that n is independent of the 
direction i. 

Whereas the correlation length ^ = C,eq{i) depends on the long-distance behaviour of 



G(r, t), this is not the case for the line-zero density n^ero- In our Gaussian approximation [^6 
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|27| of the previous section it is determined completely by the short- distance behaviour 
of G{r, t) as 



- -i/"(0.<). where /(,..t)^§M 

Some caution is necessary. Since thermal and quantum fluctuations will give rise 
to zeroes on all scales, neither G{0,t) nor its derivatives exist because of ultraviolet 
divergences. For the moment, we put in a cutoff / = 0(^o) by hand, as 

Gi{r,t)= J fke'^-''G{k,t)e-'''^\ (5.52) 

We note that the inclusion of a cut-off does not affect the long-distance correlation ^, 
which depends essentially on the weighting of the nearest singularities of G{k, t) in the 
complex fc-plane. 

In the linear regime everything is calculable. If we define the line-zero separation 

^zeroil) by 

it is apparent that ^zero{t) has little, if anything, to do with ^ directly. 

For non-Gaussian fields the situation is much more complicated. However, as long 
as there is a dominant wavenumber ko in G{k;t) this sets a length scale ^ ~ k^^ that 
characterises vortex separation. At the level of density calculations, as distinct from 
length distributions, this is essentially all that is needed. 

We can already anticipate a mismatch between the defect density and the long-range 
correlation length that could derail the Kibble-Zurek predictions. 



5.3 Annular experiments require a new correlation length 

We conclude with a brief discussion of annular experiments, that have no counterpart in 
QFT. To the extent that an annulus can be treated as a one-dimensional system there 
has been good numerical work. However, to date there are no reliable calculations for 
quenching in an annulus that take realistic boundary conditions into account. Nonethe- 
less, we have some hints as how to proceed. Consider a circular path in the bulk fluid 
(in the 1-2 plane), circumference C, the boundary of a surface 5*. For given field con- 
figurations 0a (x) the phase change 9c along the path can be expressed as the surface 
integral 

9c = 271 f (fx p3(x), (5.54) 
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Again we quench from an initial state with no rotation. The variance in the phase 
change around C, A9c is determined from 



{Mcf = 47r2 / d^x [ d'y {p,i^)p,iy)),. (5.55) 

Jx€S Jy€S 



On using the conservation of charge it is not difficult to show, from the results of |2g, ^ 
that, for Gaussian fields, with momentum cut-off at A; = A6c satisfies 



(AOcY = - f d^x I d^y CK|x-y|,t), (5.56) 



where x and y are in the plane of S, and 



C,(M)^i#f,i%L). (5.57) 



rdr\l — fi{r, t) 

Since Gi{r,t) is short-ranged Ci{r,t) is short-ranged also. With x outside S, and y 
inside S, all the contribution to (A^,^)^ comes from the vicinity of the boundary of S, 
rather than the whole area. That is, if we removed all fluid except for a strip from the 
neighbourhood of the contour C we would still have the same result. This supports 
the assertion by Zurek that the correlation length for phase variation in bulk fluid is 
also appropriate for annular flow. The purpose of the annulus (more exactly, a circular 
capillary of circumference C with radius / ^ C) is to stop this flow dissipating into the 
bulk fluid. 

More precisely, suppose that C ^ ^(t). Then, if we take the width 21 of the strip 
around the contour to be larger than the correlation length of Ci{r,t), Eq. p.56| can be 
written as 



(A9c?^-2C drr^Ci(r,t). (5.58) 
Jo 



The linear dependence on C is purely a result of Gaussian fluctuations. 

Insofar as we can identify the bulk correlation with the annular correlation, instead 
of Eq. p.lT] , we have 

The step length ^s{t) is given by 

drJ!l^0-. (5.60) 



If we quench in an annular capillary of radius / much smaller than its circumference, 
we are, essentially, coarsegraining to that scale. That is, the observed variance in the flux 
along the annulus is ttI'^Av for Av averaged on a scale /. We make the approximation 
that that is the major effect of quenching in an annulus. This cannot be wholly true, but 
it is plausible if the annulus is not too narrow for boundary effects to be important. 
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6 Nonetheless, Zurek's Predictions are Correct When 
Thermal Fluctuations are Small 



An integral part of the Zurek/Kibble predictions for Udef is their universality. The ratio 
Eq. |5.51| for n^ero gives intimations as to how this could happen, insofar as the the specific 
effects of the interactions can be subsumed into prefactors that cancel. Although, in 
the context of Gaussian fluctuations, it can only be approximate, it is part of a general 
truth that defect density only depends on very limited attributes of the fluctuation power 
spectra. 

We can see this in the models that we have introduced. The TDGL equation is highly 
dissipative, whereas the QFT equation is not, as it stands. Nonetheless, in each case the 
field will become ordered after the onset of the transition by the growth of unstable 
long wavelength modes. How this happens depends on the specifics of the equations. 
We would therefore expect the field behaviour as the transition is being implemented 
to differ in detail for condensed matter and QFT, but only weakly affecting the initial 
line-zero density, and it s abihty to describe vortices. 



6.1 TDLG condensed matter: vortex densities 



We begin with condensed matter, which we will find to be easier. As the system evolves 
away from the transition time, the free equation Eq. [4.29| ceases to be valid, to be replaced 
by the full equation 



6,(x, t) = -[-V' + e(t) + ;3|0(X, t) p]0a(x, t) + r]a{^, t), 



(6.61) 



where (3 is the rescaled coupling. 

In order to retain some analytic understanding of the way that the density is such 
an ideal quantity to make predictions for, we adopt the approximation of preserving 
Gaussian fluctuations by linearising the self-interaction as 



/),(x, t) = -[-V' + ee//(t)]0„(x, t) + r]a{x, t), 



(6.62) 



where eg// contains a (self-consistent) term O(/3(|0p)). Additive renormalisation is nec- 
essary, so that eeff ~ e, as given earlier, for t < to- 



Self-consistent linearisation is the standard approximation in non-equilibrium QFT 1 23 



but is not strictly necessary here, since numerical simulations that identify line zeroes of 
the field can be made that use the full self-interaction 11141. However, to date none ad- 



dress the questions we are posing here exactly, and until then there is virtue in analytic 
approximations provided they are not taken too seriously. 



The solution for G{r,t) is a straightforward generalisation of 4.32 



G{r,t) 



1 \3/2 



(6.63) 
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Putting in the momentum cutoff k ^ > I = lC,o = 0{C,o) of Eq. |5.52| by hand corresponds 
to damping the singularity in G{r,t) at r = as[|l5| 

making Gi{0,t) finite. We stress that, for t ~ Iq, the correlation length ^ remains O(^), 
independent of /. 

Assuming a single zero of €eff{t) at t = to, a.t r = the exponential in the integrand 
peaks at r = f = 2{t — to). Expanding about f to quadratic order gives 

For times t > e^Tq we see that, as the unfreezing occurs, long wavelength modes with 
k'^ < t/TQ — eo grow exponentially. 

The effect of the back-reaction is to stop the growth of ^^(0, t) — Gi{0, to) = (|0P)t — 
(|0p)o at its symmetry-broken value f5~^ in our dimensionless units. A necessary condi- 
tion for this is lim„_^oo ^eff{u) = 0. That is, we must choose 



(t) = e(t) + P{Gi{0, t) - Gi{0, to)), (6.66) 



thereby preserving Goldstone's theorem. 

Beyond that, what is remarkable in this approximation is that the density of line ze- 
roes uses no property of the self- mass contribution to eeff{t), self-consistent or otherwise. 
With 

-G"(0,t) = ^/ — f(t-r/2)(-^) e-/o '^^^^/i^^-^/^). (6.67) 
2 JO T xAtit J 



all pref actors in UzeTo cancel^, to give[|l^, 0] 

/,x _ 1 Jo (r+/2)5/2 e 

4vr e-(-^(*-o))V4P ^^"^^^ 

on using the definition tq = P in natural units. 

At t = to both numerator and denominator are dominated by the short wavelength 
fluctuations at small r. Even though the field is correlated over a distance ^ ^ / the 
density of line zeroes Uzero = 0{l^^) depends entirely on the scale at which we look. In 
no way would we wish to identify these line zeroes with prototype vortices. However, as 
time passes the peak of the exponential grows and n^ero becomes increasingly insensitive 
to /. How much time we have depends on the magnitude of /3, since once G'(0,t) has 
reached this value it stops growing. Since G(0,t) = 0(exp(((t - to)/t)2) at early times 

^Our ignoring prefactors in was fortuitous, leaving our conclusions obtained there unaffected. 
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the backreaction is implemented extremely rapidly. We can estimate the time t* at which 
this happens by substituting e(u) for eeffiu) in the expression for Gi{Q,t) above. 

For t > t* the equation for nzero{t) is not so simple since the estimate above, based 
on a single isolated zero of eeff{t), breaks down because of the approximate vanishing of 
eeff(t) ioT t > t*. A more careful analysis shows that Gi{0, t) can be written as 

where G{t, t) has the same peak as before at r = 2{t — to), in the vicinity of which 

but G{T,t) = 1 for r < 2(t — t*). Thus, for tq 3> tq, Gi{0,t) can be approximately 
separated as 

Gi{0, t) ^ Gf^{t) + G'^^(t), (6.71) 

where _ 

G^^{t) = f{t) / dT/\A7T{T + P)f^ (6.72) 
Jo 

describes the scale- dependent short wavelength thermal noise, proportional to tempera- 
ture, and 

describes the scale- independent, temperature independent, long wavelength fluctuations. 
A similar decomposition Gffi{0,t) ^ G//f^(t) + Gff^^{t) can be performed as 

G//\'^{t) = 27Tf{t) dT/\A7T{r + F)f^ (6.74) 
Jo 

and 

^"w= (8.(«-;„))v. L 

In particular, G//^^{t)/G^^{t) = 0{t-^). 

Firstly, suppose that, for t>t*, G^^{t) > G^^ {t) and Gtt^^it) > GttY^it), as would 
be the case for a temperature quench T{t) 0. Then, with little thermal noise, we have 
widely separated line zeroes, with density n^ero(t) ~ -Gn^^{t) /2TiG^^{t). With dn^ero/dl 
small in comparison to Uzero/^ at / = we identify such essentially non-fractal line- zeroes 
with prototype vortices, and n^ero with n^e/- Of course, we require non-Gaussianity to 
create true classical energy profiles. Nonetheless, the Halperin-Mazenko result may be 



well approximated for a while even when the fluctuations are no longer Gaussian [20 



This is supported by the observation that, once the line zeroes have straightened on 
small scales at t > t*, the Gaussian field energy, largely in field gradients, is 

F^iJ^d'x^iV^af) = -VG"iO,t), (6.76) 
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where V is the spatial volume. This matches the energy 

E ^ Vn,e/(t)(27rG'(0,t)) = -VG"{0,t) (6.77) 

possessed by a network of classical global strings with density Uzero, in the same approx- 
imation of cutting off their logarithmic tails. 
For times t > t* 



the solution to Vinen's equation [p9 



- -^'^""^^ro. (6.79) 

where X2 = 47r/ir = ATrh/aoTQ = Att if, as earlier, we motivated tq from the Gross- 
Pitaevskii equation, in which aoTo = ^- More realistically, we find X2 > 47r both for '^He 
and ^He. Taking tq ^ 8.0 x 10~^^s and ~ 5.6A) in the mean-field approximation for 
'^He gives ^2 ~ 5 x An. For ^He (with ~ 77 nm, tq ~ 1 ns) we find X2 ~ 10 x 4tt. 

This decay law is assumed in the analysis of the Lancaster experiments, in which the 
density of vortices is inferred from the intensity of the signal of scattered second sound. 
RG improvement leaves '^He unchanged, but for '^He it redefines X2 to X2i^ — T/Tc)~^^^ > 

X2- 

The first problem for the Lancaster experiments is that this makes an already large X2 
even larger. In the attenuation of second sound the signal to noise ratio is approximately 
0(l/x2^)- The empirical value of X2 used in the Lancaster experiments is not taken 
from quenches, but turbulent flow experiments. It is suggested [p!^] that X2 ~ 0.005, a 
good three orders of magnitude smaller than our prediction above. Although the TDLG 
theory is not very reliable for ^ife, if our estimate is sensible it does imply that vortices 
produced in a temperature quench decay much faster than those produced in turbulence. 
^He experiments provide no check. 

This is only one of our worries. We shall argue that, for early time at least, thermal 
fluctuations are large in the Lancaster experiments. However, for ^He, with negligable 
UV contributions, we estimate the primordial density of proto-vortices as 



"^^eroit*)^———--., —, (6.80) 




in accord with the original prediction of Zurek. Because of the rapid growth of G{0,t), 
{t* — to) /i = p > 1 = 0(1). With p behaving as (ln(l//?))^/^ there is very little variation. 
For ^He quenches p ^ 5 (and for '^He quenches p ~ 3). We note that the factor^] of 
P = 87rp gives a value of / = 0(10), in agreement with the empirical results of |jTO| and 
the numerical results of 



^An errant factor of 3 appeared in the result of [Q 

^The temperature quench of the latter is somewhat different from that considered here, but should 



still give the same results in this case 
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6.2 Path integrals 

Finally, we note that G{r, t), and hence its derivatives, can be expressed as path integrals 
for the diffusion of a particle in a time-dependent potential eg//. Thus, for example. 



^/;/ N T dr r ^ \ /"^ , i.dx,. 



1 ,(ix. 



(6.81) 



In Eg . [6 .811 the summation is over closed paths for the particle (mass 2, in our units) 



traversed in time r. We note that, if e = eo is fixed, then 

T dr 



G"(0,t) = - \ —i Pxe-^(")^o (6.82) 

2 JO T J 



where L{j) = Jq dr' Jii^ij') is the length of the loop. Eq. |6.82 is just the free energy of 



a 'free gas' of fluctuating relativistic Brownian loops of varying lengths and shapes with 
tension eo. In conventional statistical field theory the vanishing of eo in Eq. |6.82| signals 
a phase transition due to the proliferation of loops. The situation here is somewhat 
different in that the loops have a tension that varies with both the external time t and 
r. The onset of the transition is signaled by this tension vanishing at some point on the 
loops. 

7 For Slow Pressure Quenches or a Large Ginzburg 
Regime Thermal Fluctuations Cannot be Ignored 

7.1 ^He experiments 

The situation for the Lancaster '^He experiments is complex, since they are pressure 
quenches for which the temperature T is almost constant at T ^ Tc. Unlike temperature 
quenches fl^, |3T[| , thermal fluctuations here remain at full strength]^. The necessary time- 



mdependence of G^^{t) for t > t* is achieved by taking e^ffiu) = 0{u~^). In consequence, 
as t increases beyond t* the relative magnitude of the UV and IR contributions to Gz(0, t) 
remains approximately constant. Further, since for t = t*, 

this ratio is the ratio ai t = t*. 

Nonetheless, as long as the UV fluctuations are insignificant at t = t* the density of 
line zeroes will remain largely independent of scale. This follows if G//^^{t*) ^ G//y^{t*), 
since Gffi{0,t) becomes scale-independent later than Gi{0,t). In we showed that this 
is true provided 

(Vj-o)(l - Tg/T,) < Gtt\ (7.84) 

^Even for ^He, T/Tc never gets very small, and henceforth we take T = Tc in Gi{0,t) above 
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where C = 0(1) and Tc is the Ginzburg temperature. With tq/tq = O(IO^) and 
(1 — Tg/Tc) = O(10~^^) this inequahty is well satisfied for a linearised TDLG theory 
for '^He derived^ from the full TDGL theory[^, but there is no way that it can be 
satisfied for ^ife, when subjected to a slow mechanical quench, as in the Lancaster 
experiment, for which tq/tq = O(10^°), since the Ginzburg regime is so large that (1 — 
Tg/Tc) = 0(1). Effectively, the ^ife quench is nineteen orders of magnitude slower 
than its ^He counterpart. It is satisfying to find the Ginzburg regime reappear as an 
essential ingredient in undermining universality. In particular, as in the original proposal 
of Kibblefll], large thermal noise inhibits vortex production. 

When the inequality is badly violated, as with ^ife for slow pressure quenches, then 
the density of zeroes n^ero = 0(/~^) after t* depends exactly on the scale I at which we 
look and they are not candidates for vortices. Since the whole of the quench takes place 
within the Ginzburg regime this is not implausible. However, it is possible that, even 
though the thermal noise never switches off, there is no more than a postponement of 
vortex production, since our approximations must break down at some stage. The best 
outcome is to assume that the effect of the thermal fluctuations on fractal behaviour is 
diminished, only leading to a delay in the time at which vortices flnally appear. Even 
if we suppose that n^ero above is a starting point for calculating the density at later 
times, albeit with a different to, thereby preserving Vinen's law, we then have the earlier 
problem of the large X2 = 0{p). In the absence of any mechanism to reduce its value 
drastically, this would make it impossible to see vortices. As a separate observation, we 
note that the large value of p in the prefactor of n^eTo is, in itself, almost enough to 
make it impossible to see vortices in ^ife experiments, should they be present. This will 
be pursued elsewhere. 

In summary, this work suggests that for slow pressure quenches in ^ife, we see no 
well-deflned vortices at early times because of thermal fluctuations, and it is plausible 
that, if we do see them at later times, there are less than we would have expected because 
of their rapid decay and their initial low density. The situation is different for ^He. That 
the density of such vortices as appear may agree with the Zurek prediction is essentially 
a consequence of dimensional analysis, given that the main effects of the self-interaction 
have a tendency to cancel in the counting of vortices, without introducing new scales. 
However, a numerical simulation that goes beyond the Gaussian approximation that is 
speciflcally tailored to the Lancaster parameters is crucial if we are to understand the 
late-time behaviour and see if this suggestion can be sustained. We shall pursue this 
elsewhere. 

7.2 TDLG: annular physics 

We saw for vortex formation in '^He that the dependence of the density on scale makes 
the interpretation of observations problematic. This is not the same if the ^He is conflned 
to an annulus, since the annulus itself provides the coarse-graining. That the incoherent 

^"Ignoring the position-dependent temperature of [|| 
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depends on its radius / is immaterial. The end result is that 




The time t* for ^ife at which we evaluate v, when = ao/P, depends weakly on /, 

varying from about 3t to At as / varies from I = <^ ^ to I = 10^. For / > 4^ the scale 
at which the coarse-grained field begins to occupy the ground states becomes largely 
irrelevant. 

On expanding /'^(r, t; /)/(! — /^(r, t; /)) in powers of r and keeping only the forward 



peakp^ 



1 4G2^3G3 G2 \-V2 

' ' , (7.86) 



^,{t;l) 9GiV20G'2 12^1, 

where Gn{r,t) is the {2n — 2)th derivative of G with respect to r at r = (and hence 
proportional to the 2?T,th moment of G{k,t)). 

If, as suggested by Zurek, we take Z = O(^) we recover Eq. p.lT] qualitatively, although 
a wider bore would give a correspondingly smaller flow. 



We assume 32 that 



2/>&//(M) = (|;|-^)"''. (7.87) 

since otherwise the correlations in the bulk fluid from which we want to extract annular 
behaviour are of longer range than the annulus thickness. The effect of this is to reduce 
the flow velocity for narrower annuli. The effect is largest for small radii / < ^, for which 
the approximation of trying to read the behaviour of annular flow from bulk behaviour 
is most suspect. 

Once / is very large, so that the power in the fluctuations is distributed strongly 
across all wavelengths we recover our earlier result, that C,s{t*]l) = 0{l). However, the 
change is sufficiently slow that annuli, significantly wider than ^, for which experiments 
are more accessible, will give almost the same flow as narrower annuli. This would 
seem to extend the original Zurek prediction of Eq. |3.17| to thicker annuli, despite our 



expectations for incoherent flow. However, we stress again that caution is necessary, 
since in the approximation to characterise an annulus by a coarse-grained ring without 
boundaries we have ignored effects in the direction perpendicular to the annulus. In 
particular, the circular cross-section of the tube has not been taken into account. One 
consequence of this is that infinite (non-self-intersecting) vortices in the bulk fluid have 
no counterpart in an annulus. 

Since Av only depends on ^7^^^ it is not sensitive to choice of / > 2^ at the relevant t. 



Given all these approximations our final estimate|32| is (in the cm/sec units of Zurek 



Av ^ 0.2{rQ[fis])-''/y^C[cm] (7.88) 



28 



for radii of 2^ — 4,^, tq of the order of milliseconds and C of the order of centimetres. 
1/ = 1/2 is the mean- field critical exponent above. In principle v should be renormalised to 
1/ = 2/3, but the difference to Aw is sufficiently small that we shall not bother. Given the 
uncertainties in its derivation the result Eq. |7.88| is indistinguishable from Zurek's0] (with 
prefactor 0.4), but for the possibility of using somewhat larger annuli. The agreement 
is, ultimately, one of dimensional analysis, but the coefficient of unity could not have 
been anticipated easily, given the small prefactors of Eq.6.20. How experiments can be 
performed, even with wider annuli, is another matter. 



8 QFT: The Appearance of Structure 

It is because the formation of defects is an early-time occurrence that it is, in part, 
amenable to analytic solution. Again we revert to the mode decomposition of Eq. [4.42| . 



The field becomes ordered, as before, because of the exponential growth of long- wavelength 
modes, which stop growing once the field has sampled the groundstates. What matters is 
the relative weight of these modes (the 'Bragg' peak) to the fiuctuating short wavelength 
modes, since the contribution of these latter is very sensitive to the cutoff I at which we 
look for defects. Only if their contribution to Eq. |2.14| is small when field growth stops 



can a network of vortices be well-defined at early times, let alone have the predicted 
density. Since the peak is non-perturbatively large this requires small coupling, which 
we assume. However, there is a problem in that, in the absence of explicit damping of 
the type seen in FRW universes, rescattering of modes can rapidly undo the early-time 
appearance of structure. Thus, while we take small coupling, we do not take very small 
couplings of the magnitude (e.g. 10^^^) associated with infiationary models. 



8.1 Mode Growth v Fluctuations: The free roll (continued) 



We begin by extending the analysis of Section 4.2 to later times, still in the approximation 
of a free roll. This needs care for slow quenches since the backreaction serves to hold the 



field in the vicinity of the intermediate groundstates |0 | 

-m2(t) M^At 



00 (i^) where, now 



\it) 



X 



59) 



where At = t — to = t — coTq as before. Nonetheless, the free roll provides a basis for the 
more realistic picture. 

Prior to the completion of the quench at At = tq, the mode equation ( 4.42 ), now of 
the form 

;iT^ + k' X^W = 0, (8.90) 

is exactly solvable, as we saw earlier. 
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For this section it is convenient to redefine the origin of time at t = to? whereby we 
can drop the prefix A. 

We are primarily interested in the exponentially growing modes that appear when 

nl{t) = -k' + ^>0. (8.91) 

For fixed k this occurs when t > tj. = Tqk"^ /M"^. 

The WKB solution is adequate for our purposes. The coarse-grained Gi{r;t) can be 
written as Gi{r; t) = G'^^P{r; t) + Gf'^'^{r; t) where 

i\k\<kt \kt 

has exponentially growing long wavelength modes and 



G^^nr;t) ^ ^ / #fe .f e^''-\ath„{Su{t)) + a-J.,/,{S,m' (8-92) 



has short wavelength oscillatory modes. In both cases 

2 M 



S,{t) = f dt'\Qk{t')\ = -^\t-t,\'/'. (8.94) 
Jt- i ^rrn 



As in the case of condensed matter previously, we have coarse-grained the field by 
introducing a simple cut-off at A; = A = 0(M), or / = A~^. For fixed t the dividing 
monentum is k\ = M'^I/tq. Provided we are far from the transition we have incorporated 
the initial data into the in ( p.92| ) and ( p. 931) . The normalisation factor Tq/M'^ has 
been made visible. The remaining have no A dependence. Since G°^^{0]t) = 0{ToM) 
for — t = 0{tq) the have no tq dependence. 

For large t the integrand in ( p.92| ) will be peaked at some A;o(t) — as t — *• cxo, once 
the angular integrals have been performed. Assuming that /co(^) ^ kf the upper bound 
in the integral can be dropped and \k^ — approximated by \k^\, knowing that there 
is no singularity at = kf. With nothing to stop + behaving like a nonzero 
constant in the vicinity of = 0, it can be treated as slowly varying and the integral 
approximated as [p3 



^ ' ' M^\t ) J\k\<M ^ 

oc ( '^"e^Mtyys^^-rVeit) (8,95) 

M\m{t)\\^J ^ ' 

on performing the k"^ expansion of the exponent, where 

e{t) = (8.96) 

Although, like Eq. p.95| , the expression Eq. p.96| is not supposed to be valid for small t, it 
does embody the Kibble freezeout condition Eq. |2.10| in satisfying C,{t) = 0(1). 
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8.2 Comparison with Kibble's results: First guess 

The calculations we did on CM systems showed that, although the freeze- in time i was 
not relevant for the appearance of stable vortices, their subsequent density is determined 
by the scales at this time in the absence of strong fluctuations. We shall see that the 
same could be true here. 

The calculations above were for a free roll. Let us suppose, provisionally, that the 
backreaction exerts its influence over such a short time that, in effect, it is if it were 
an instantaneous brake to domain growth. The provisional freeze-in time t* is then, 
effectively, the time it takes to reach the transient groundstate (pl{t) = — m^(t)/A. That 
is, G(O;r)=O(02(r)), giving 

Neglecting Inln terms gives 

2/3 



Mt* 

2 



-(MrQ)i/3(^ln(l/A) + ln(MrQ)j . (8.98) 



The second term in the bracket of (|8.98|) can be ignored to a good approximation provided 
Mtq < (1/A), preferably by a large margin, and we assume that this is so. In terms of 
the causal time t the relation is just 

Mt* ~ Mt(ln(l/A))2/l (8.99) 

That is, the freeze in time t* is (qualitatively) larger than the causal time t. As far as the 
separation of scales is concerned, we have the same effect qualitatively if we had taken 
t* as the time for the field to reach the final ground state as |0^| = M^/A, rather than 
the provisional ground states 0o(^)- Nonetheless, we shall remain with (|8.97|) . 

Whatever, the necessary condition that t* < tq (since the previous calculations are 
for t < Tq) reduces to 

Mtq > ln(l/A). (8.100) 

where, and henceforth, we neglect coefficients of 0(1). That is, the quench time should 
be longer than the freeze-in time for the instantaneous quench, the time it takes the field 
to sample the ground-state in a free-roll^, but shorter than the equilibriation time (or 
powers of it that preserve the separation of scales). 

At this qualitative level the correlation length at the spinodal time is 

M^^\t*) ^ (MrQ)2/3(ln(l/A))^/l (8.101) 

The effect of the other modes is larger than for the instantaneous quench, giving, at 
t = t* 
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In general, we recover the results of the instantaneous quench if we set Mtq — 0(ln(l/A)) 
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The error term E = 0(A^/^(MrQ)^/^(ln(l/A)) ^^^) is due to oscillatory modes, sensitive 
to the cut-off. In mimicry of Eg . p .141 it is helpful to rewrite Eq. |8.102] as 







Mi ' 





(ln(l/A))-^/^[l + E]. 



(8.103) 



in terms of the scales tq = = M~'^. The first term in Eq. |8.103] is the Kibble estimate 
of Eq. |2.14| , the second is the multiplying factor, rather like that in Eq. |6.80| , that yet 
again shows that estimate can be correct, but for completely different reasons. As for 
condensed matter, the dependence on the interaction strength is only through a power 
of the logarithm of (1/A))^^''^. 

The third term shows when it can be correct, since E is also a measure of the sensi- 
tivity of nzero to the scale at which it is measured. The condition i?^ <^ 1, necessary for 
a vortex network to be defined, is then guaranteed]^ if 

(8.104) 



{tq/to)\1-Tg/T,)<C, 



where C = 0(1), on using the relation (1 — Tg/Tc) = 0(A). This is the QFT counterpart 
to Eq. |7.84| and concurs again, in principle, with Kibble's earlier argument that large 
thermal fluctuations inhibit the appearance of vortices. 

The easiest way to enforce E <^ 1 and Mtq > ln(l/A) is to take Mtq = ln(l/A)", 
for a > 1. The effect in Eq. 8.103| is merely to renormalise the critical index. Of course, 
the Kibble prediction Eq. |2.14| was only an estimate. Although it is good qualitatively 
it is misleading when considering definition of the network since a simple calculations 
shows that, at time t, the string density is still totally sensitive to the definition of 
coarse-graining. 

Finally, suppose that this approach is relevant to the local strings of a strong Type-II 
U{1) theory for the early universe, in which the time-temperature relationship tT"^ = 
TMpi is valid, where we take T = 0(10"^) in the GUT era. If G is Newton's constant 
and /i the classical string tension then, following 0, Mtq ~ 10~^A^/^(G/i)~^''^. The 
dimensionless quantity Gfi ~ 10~^ — 10~^ is the small parameter of cosmic string theory. 
A value A ~ 10^^ gives Mtq ~ (Mt*)",a ~ 2, once factors of tt, etc. are taken into 
account, rather than Mtq ~ 1/A, and the density of Eq. |8.103] may be relevant. 



8.3 Backreaction in QFT 

To improve upon the free-roll result more honestly, but retain the Gaussian approximation 
for the field correlation functions, the best we can do is adopt a mean-field approximation 
along the lines of p3|, p5[, as we did for the CM systems earlier. As there, it does have 



the correct behaviour of stopping domain growth as the field spreads to the potential 
minima. As before, only the large-iV expansion preserves Goldstone's theorem. 

G{r; t) still has the mode decomposition of ( [4.44|) , but the modes xt ^o"^ satisfy the 
equation 

^ + k2 + m2(t) + A($2(0))t 



(8.105) 
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where we have taken N = 2. Because A(/)^ theory is not asymptotically free, particularly 
in the Hartree approximation, the renormalised A coupling shows a Landau ghost. This 
means that the theory can only be taken as a low energy effective theory. 
The end result is||23|, on making a single subtraction at t = 0, is 



^ + + m\t) + A J fpC{p)[x;{t)x,{t) - 1] 



which we write as 



XfeW = 0. 



(8.106) 



^.107) 



On keeping just the unstable modes in ($^(0))t then, as it grows, its contribution to 
( p.l06| ) weakens the instabilities, so that only longer wavelengths become unstable. At t* 
the instabilities shut off, by definition, and oscillatory behaviour ensues. Since the mode 
with wavenumber > stops growing at time t\ < t*, where /i^(t^) = k^, the free- roll 
density at t* must be an overestimate. 

An approximation that improves upon the WKB approximation [HI] is 



Xk{t) 



TtM \1/2 



exp 



dtn{t) 



(8.108) 



,2fife(r/). 

when rj = M(t^ — t) > is large, and ^l{t) = fi'^{t) — k^. On expanding the exponent in 
powers of k and retaining only the quadratic terms we recover the WKB approximation 
when /i(t) is non-zero. 

The result is that the effect of the back- reaction is to give a time-delay At to t*, 
corresponding to a decrease in the value ko{t) at which the power peaks of order 

The backreaction has little effect for times t < t*. For t > t* oscillatory modes take over 
the correlation function and we expect oscillations in G{k;t). 

In practice the backreaction rapidly forces fi'^{t) towards zero if the coupling is not too 



small For couplings that are not too weak the end result is a new power spectrum, 
obtained by superimposing oscillatory behaviour onto the old spectrum. As a gross 
oversimplification, the contribution from the earlier exponential modes alone can only be 
to contribute terms something like 

T 



G{r;t) oc 



^4M{t*)3/V3 



X 



M\m{t*)\ 
cosk{t - t*) + 

-1/2 



'ke 



ik.x 



|k|<Af 

n{k) - w'jk) 
k 



sin k{t-t*] 



(8.110) 



to G, where = M{t* - tfc)^/V^Q and W = l/4(t* - tk). Fortunately, the details are 
almost irrelevant, since the density of line zeroes is independent of the normalisation, 
and only weakly dependent on the power spectrum. 
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The k = mode of Eq. |8rTT0| encodes the simple solution xa:=o(^) = a+bt when /x^ = 0. 
As observed 0] by Boyanovsky et al. this has built into it the basic causality discussed 
by Kibblep. Specifically, for r,t ^ oo, but r/t constant (7^ 2), 

G(r,t) ^ -e(2t/r- 1). (8.111) 
r 

It has to be said that this approximation should not be taken very seriously for large t, 
since we would expect rescattering to take place at times At = 0(1/A) in a way that the 
Gaussian approximation precludes. 

Whatever, it follows directly that this causality, engendered by the Goldstone particles 
of the self-consistent theory, has little effect on the density of line- zeroes that we expect 
to mature into fully classical vortices. If Eq. ^.l04] is not satisfied, it is difficult to imagine 



how clean vortices, or proto- vortices, can appear later without some additional ingredient. 
8.4 Really slow quenches 

Finally, consider slowing the quench until the WKB approximation manifestly breaks 
down. This is almost certainly the case if the growing modes catch up with the moving 
minima within the Ginzburg regime by time 

tc = 0{\tq). (8.112) 

If we take ( ^.97| ) seriously it follows that domain growth has stopped by time to provided 
Mtq is a power of A~^. Equation (|8.97|) only makes sense if Mtc 3> 1 i.e Mtq 3> A~^. As 



an example, let us take Mtq = A~'^/^. In this case the dynamical correlation length ^(tc) 
of ( p.96| ) is 0(1/MA^/^), equal to the equilibrium correlation length ^eqitc) = |"^~^(^g)|5 
suggesting the string density no = 0{M'^/X) that follows from ( p.6|) , very much smaller 



than that of ( p.l03| ) above. However, a simple check shows that, in this case, 

I driiitc) 



riiitc) dl 



l=M~ 



Oil). (8.113) 



Despite the low density of line zeroes at the Ginzsburg temperature, they do not provide 
a stable network of proto-vortices. 



9 Conclusions 

We examined the Kibble /Zurek predictions for the onset of phase transitions and the 
appearance of defects (in particular, vortices or global cosmic strings) as a signal of the 
symmetry breaking. Our results are in agreement with their prediction Eq. |2.13| as to the 
magnitude of the correlation length at the time the transition truly begins, equally true 
for condensed matter and QFT. 
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However, this is not simply a measure of the separation of defects at the time of their 
appearance. The time i is too early for the field to have found the true groundstates 
of the theory. We believe that time, essentially the spinodal time, is the time at which 
proto-vortices can appear, which can later evolve into the standard classical vortices of 
the theory. 

Even then, they may be frustrated by thermal field fluctuations. In TDLG condensed 
matter thermal noise is proportional to temperature. If temperature is fixed, but not 
otherwise, as in the pressure quenches of '^He this noise can inhibit the production of 
vortices, although there are other factors to be taken into account (such as their decay 
rate). On quenching from a high temperature in QFT there are always thermal fluctua- 
tions, and these can also disturb the appearance of vortices. The condition that thermal 
fluctuations are ignorable at the time that the field has achieved the true ground-states 
can be written 

irQ/Toni-TG/T,)<C, (9.114) 

where 7 = 1 for condensed matter and 7 = 2 for QFT. C = 0(1). 

This restores the role of the Ginzburg temperature Tq that the simple causal argu- 
ments overlooked. As in the earlier arguments, large fluctuations inhibit vortex produc- 
tion, but in their absence the equilibrium correlation length is not the relevant quantity. 
Quenches in "^He provide the major example for which Eq. |9.114] is not satisfied. However, 



when Eq. |9.114" is satisfied, the Kibble/Zurek prediction is recovered for the density of 



line zeroes, now seen as potential vortices. That no new scales (up to logarithms) are 
introduced by the interaction is a reflection of the fact that line-zero density only uses 
limited properties of the power spectrum of fluctuations (in the Gaussian approximation 
at least). 

What happens at late time is unclear, although for TDLG numerical simulations can 
be performed (but have yet to address this problem exactly). On the other hand, not only 
is the case of a single self-interacting quantum scalar field in flat space-time a caricature of 
the early universe, but it is extremely difficult to go beyond the Gaussian approximation. 
To do better requires that we do differently. There are several possible approaches. One 
step is to take the FRW metric of the early universe seriously, whereby the dissipation 
due to the expansion of the universe can change the situation dramatically ||36|| . Other 
approaches are more explicit in their attempts to trigger decoherence explicitly, as we 
mentioned earlier. Most simply, one treats the short wavelength parts of the field as 
an environment to be integrated over, to give a coarse-grained theory of long- wavelength 
modes acting classically in the presence of noise. However, such noise is more complicated 
than in TDLG theory, being multiplicative as well as additive, and coloured [p^, p2| , pT| . 
This is an area under active consideration, but we shall stop here. 

I would like to thank Alasdair Gill, Tom Kibble, Wojciech Zurek and Grisha Volovik 
for fruitful discussions. There is much work in this area and we have not included all 
references. We apologise to any authors who we may have missed inadvertently. This 
work is the result of a network supported by the European Science Foundation. 
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